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15 ABSTRACT: In order to compare stand structure and woody species diversity of subtropical 

16 evergreen broadleaf forests along a latitudinal thermal gradient of the Ryukyu Archipelago, 

17 tree censuses in a 750 m 2 plot in Okinawa Island and a 400 m 2 plot in Ishigaki Island were 

1 8 performed. The number of layers increased along a latitudinal thermal gradient from four in 

19 the forest of Okinawa Island to five in the forest of Ishigaki Island. The values of Shannon's 

20 index H' and Pielou's index /' tended to increase from the top layer downward in the forest 

21 of Okinawa Island, whereas in the forest of Ishigaki Island, these values tended to increase 

22 from the bottom layer upward. High woody species diversity depended on small-sized trees in 

23 the Okinawa forest, whereas it depended on large-sized trees in the Ishigaki forest. The 

24 woody species diversity is higher in the Okinawa forest (H'= 4.83 bit) than in the Ishigaki 

25 forest (H'= 4.36 bit). According to successively decreasing height of layers from the top 
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1 downward, the value of //'increased continuously from the top layer downward in the 

2 Okinawa forest. This increasing trend was different from the Ishigaki forest, where the value 

3 of H ' increased up to the second layer and then decreased downward. In the Okinawa forest, 

4 the expected number of species increased continuously from the top toward the bottom layer, 

5 i.e. the bottom layer contained the highest potential number of species (65). However, in the 

6 Ishigaki forest, it increased from the top to the fourth layer and then decreased to the bottom 

7 layer, i.e. the fourth layer contained the highest potential number of species (90). The floristic 

8 composition in the Okinawa forest was different from that in the Ishigaki forest in terms of 

9 similarity index, though approximately half of the species were common between them. The 

10 highest degree of similarity in floristic composition was between the second and third layers 

11 in the Okinawa forest, whereas it was between the third and bottom layers in the Ishigaki 

12 forest. The degree of similarity in floristic composition among layers was higher in the 

1 3 Okinawa forest than in the Ishigaki forest. Except the top and the bottom layer respectively in 

14 the forests of Okinawa Island and Ishigaki Island, the spatial distribution of trees was random 

15 in each layer. The degree of overlapping in the spatial distribution of trees among layers in 

16 these two forests suggested that trees in the upper two layers in the Ishigaki forest can catch 

17 sufficient light, while light can not penetrate easily to the lower three layers in both of the 

18 forests. As a result, almost species in the lower layers might be shade-tolerant in both of the 

19 forests. For both of the forests, mean tree weight of each layer decreased from the top 

20 downward, whereas the corresponding tree density increased from the top downward. This 

21 trend resembled the mean weight-density trajectory of self-thinning plant populations. 
22 

23 

24 KEY WORDS'. Architectural stratification; Floristic similarity; Ryukyu Archipelago; 

25 Subtropical evergreen broadleaf forest; Woody species diversity. 
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1 

2 

3 INTRODUCTION 

4 

5 The coastal areas of the Western Pacific from subarctic eastern Siberia to equatorial SE Asia 

6 have forest climates with sufficient rainfall, which develops a sequence of five forest 

7 formations: subarctic evergreen conifer forests, cool-temperate deciduous broadleaf forests, 

8 warm-temperate lucidophyll forests, subtropical forests and tropical rain forests (Kira, 1991). 

9 Within the Western Pacific sequence of thermal vegetation zones, the subtropical zone whose 

10 major part is covered by dry area, only a small part including a chain of islands from 

1 1 Okinawa to Taiwan and South China, is sufficiently moist to allow the development of 

12 subtropical forests. Therefore, the subtropical forests in the Ryukyu Archipelago are precious 

1 3 from a phytogeographical viewpoint. 
14 

15 A well -developed evergreen broadleaf forest exists in the northern part of Okinawa Island 

16 and in the central part of Ishigaki Island of the Ryukyu Archipelago. It is ecologically 

17 important to know how floristic composition, spatial distribution of trees, multi-layering 

1 8 structure and woody species diversity among layers change along a latitudinal thermal 

19 gradient of the Ryukyu Archipelago. 
20 

2 1 The diversity of tree species is fundamental to total forest biodiversity, because trees provide 

22 resources and habitats for almost all other forest species (Hall & Swaine, 1976; Huston, 1994; 

23 Whitmore, 1998; Huang et al., 2003). Measures of species diversity play a central role in 

24 ecology and conservation biology. The most commonly employed measures of species 

25 diversity are the Shannon function, species richness (number of species), and evenness (the 
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1 distribution of abundance among the species, sometimes known as equitability). In addition, 

2 the spatial distribution of trees has been a major source of interest for plant ecologists because 

3 of its potential role in explaining the coexistence of tree species in species-rich forests 

4 (Bunyavejchewin et al., 2003). 
5 

6 The degree of canopy multi-layering and the woody species diversity increase along a 

7 latitudinal thermal gradient from higher latitudes to the tropics (Hozumi, 1975; Yamakura, 

8 1987; Kira, 1991; Ohsawa, 1995; Kimmins, 2004). The canopy multi-layering structure, i.e. 

9 architectural stratification, is an important factor in maintaining higher woody species 

10 diversity (Roberts & Gilliam, 1995; Lindgren & Sullivan, 2001). However, there is a dearth of 

1 1 studies reporting the effect of the architectural stratification on floristic composition, woody 

12 species diversity and spatial distribution of trees in the subtropical evergreen broadleaf forest 

13 in Okinawa Island. It is evident that there have apparently been no studies from such a point 

14 of view in the subtropical evergreen broadleaf forest in Ishigaki Island. The purposes of this 

15 study were, therefore, to distinguish architectural stratification, to quantify woody species 

16 diversity, floristic composition and spatial distributions of trees on the basis of the 

17 architectural stratification and to compare the above parameters between the forests in the 

1 8 northern part (Okinawa Island) and the southern part (Ishigaki Island) of the Ryukyu 

19 Archipelago. 
20 

21 

22 MATERIALS AND METHODS 

23 

24 Study sites and sampling 

25 



4 



Downloaded from http://biorxiv.org/on September 18, 2014 

1 Two study sites were respectively selected in Okinawa Island and Ishigaki Island of the 

2 Ryukyus Archipelago. One is located in a subtropical evergreen broadleaf forest at Mt. 

3 Yonaha (26°45 1ST and 128°10Ti), the northern part of Okinawa Island. The bedrock is 

4 composed of silicate having soil pH of 4.35 (Alhamd et al., 2004; Feroz et al., 2006). The 

5 other is located in a subtropical evergreen broadleaf forest at Mt. Omoto (24°25'03"N and 

6 124°11 '17"E), the central part of Ishigaki Island. The bedrock is composed of silicate having 

7 soil pH of 4.55. 
8 

9 For the former site, a sampling plot of area 750 m 2 (25 m x 30 m) was established and divided 

10 into 120 quadrats (2.5 m x 2.5 m). The slope, altitude and direction of the plot were 24.5°, 250 

11 m above sea level and NWW, respectively. For the latter site, a sampling plot of area 400 m 2 

12 (20 m x 20 m) was established and divided into 64 quadrats (2.5 m x 2.5 m). The slope and 

1 3 altitude and direction of the plot were 2 1 .4°, 200 m above sea level and SW, respectively. 

14 Woody plants over 0. 10 m in height for the former site and all woody plants for the latter site 

15 were numbered. They were identified to species according to the nomenclature of Walker 

16 (1976). Tree height H (m) and stem diameter at a height of Hl\0 D 0 .ih (cm) were measured. 
17 

18 

19 Climate 

20 

21 The climatic data for the year of 2003-2007 were collected from Nago Meteorological Station 

22 for the northern part of Okinawa Island and Maezato Meteorological Observatory for the 

23 central part of Ishigaki Island. Mean monthly minimum temperature and mean monthly 

24 maximum temperature are respectively 16.2 ± 0.5 (SE) °C in January and 29. 1 ± 0.2 (SE) °C 

25 in July in Okinawa Island, and 18.2 ± 0.4 (SE) °C in January and 29.2 ±0.3 (SE) °C in July in 

26 Ishigaki Island. Mean annual temperature is 22.9 ± 0.3 (SE) °C in Okinawa Island and 23.9 ± 
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1 0.3 (SE) °C in Ishigaki Island. The warmth index is 214.2 + 0.5 (SE) °C month in Okinawa 

2 Island and 227.2 ±1.3 (SE) °C month in Ishigaki Island, which values are within the range of 

3 1 80 to 240 °C month of the subtropical region defined by Kira (1977). Mean monthly rainfall 

4 in both of the islands is over 100 mm throughout the year except for 84 ± 22 (SE) mm in 

5 February in Okinawa Island and 93 ± 19 (SE) mm in December in Ishigaki Island. Mean 

6 annual rainfall is 2050 ± 182 (SE) mm yr" 1 in Okinawa Island and 1942 + 159 (SE) mm yr" 1 in 

7 Ishigaki Island. Typhoons with strong winds and rains frequently strike the islands between 

8 July and October. 
9 

10 

1 1 Architectural stratification 
12 

13 The M-w diagram proposed by Hozumi (1975) was used to identify the multi-layering 

14 structure of the forest stands. Tree weight w was assumed to be proportional to Do.m 2 H and it 

15 was arranged in descending order. Average tree weight M„ from the maximum tree weight w\ 

16 to the nth tree weight was calculated using the form of M n = ^ w ; In (n= I,...., total 

i=i / 

17 number of trees AO- If the M-w diagram is constructed by plotting the values of M against the 

18 corresponding values of w on logarithmic coordinates, then some segments on the M-w 

19 diagram are formed. Each segment is related to the layer with the specific characteristics of 

20 the beta-type distribution designated by Hozumi (1971, 1975). Hozumi (1975) pointed out 

21 that the segments on the M-w diagram can be written by either of the following equations: 

22 M = Aw + B (1) 

23 M=Cw b (2) 

24 where A, B, C and b are coefficients. These functions reflect some aspect of the manners of 

25 packing trees into the three-dimensional space as realized by a forest stand. 

6 
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1 A boundary between layers was distinguished by the relationships of the first 

2 derivative Si and the second derivative 52 to tree weight w on the M-w diagram. The Si and 

3 52 are defined as follows: 

_ dlnM AlnM 

4 5,= = (3) 



5,= 



d In w A In w 
d (d\nM\ dS { 



d\nw\ d\nw 



AS, 



d\nw A In w 

6 When the moving value of 5i is plotted against the corresponding value of w, points of 

7 inflection probably appear instead of the points of intersection between the segments on the 

8 M-w diagram. Then, the second derivative 52 has extrema, whose x-coordinates can be 

9 regarded as the points of boundary between layers. 
10 

11 

1 2 Species dominance 
13 

14 Dominance of a species was defined by importance value IV (%) of the species: 



(4) 



15 IV = 



100 + 100 + -/^x 100 



X/. 



(5) 



16 where n i is the number of individuals of the z'th species, a i is the basal area at a height of 

17 ///10 of the z'th species, f i is the number of quadrats in which the z'th species appeared and Q 

18 is the total number of quadrats. 
19 

20 

21 Species-area relationship 

22 
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1 The expected number of species S c/ appeared within the number of quadrats q selected at 

2 random from the total number of quadrats Q was calculated from the equation proposed by 

3 Shinozaki (1963) (cf. Hurlbert, 1971): 



4 * ? = Z 

i=i 



1 



1 



(6) 



5 where q i is the number of quadrats in which the z'th species occurred and S is the total number 

6 of species. The S q -values were obtained for ^-values of 1, 2, 4, 8, 16, 32, 64 and 120 for the 

7 forest in Okinawa Island, and 1, 2, 4, 8, 16, 32 and 64 for the forest in Ishigaki Island. 
8 

9 

10 Floristic similarity 
11 

12 The similarity of floristic composition between layers was calculated using the following 

13 index C n (Horn, 1966; cf. Morishita, 1959; Kimoto, 1967): 

S 

14 Cn = (n A +n B )/v A /v B ( ? ) 



15 (n A = $> A , 2 /V , n B = $> B , 2 //v B 2 ) 

M / i'=l 



16 where S is the total number of species, um and zzb; are the number of individuals of the z'th 

17 species respectively belonging to Layer A and Layer B. The value of C n is 1.0 when the 

1 8 number of individuals belonging to a species is the same between the two layers for all 

19 species, i.e., floristic composition is completely the same between the layers, and is 0.0 when 

20 no common species is found between them. 

21 Equation (7) was also applied for measuring the degree of similarity in floristic 

22 composition between the Okinawa and Ishigaki forests. In this case, S is the total number of 
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1 species in the two forests, tiAi and nw are the number of individuals of the z'th species 

2 respectively belonging to the Okinawa and Ishigaki forests. 
3 

4 

5 Species diversity 

6 

7 The following two indices of Shannon's index (Mac Arthur & MacArthur, 1961) H ' (bit) and 

8 Pielou's (1969) index J' were used to measure woody species diversity or equitability 

9 (evenness): 

10 H '=i^2- (8) 

11 J '=-^- (#: ax = iog 2 s) (9) 

max 

1 2 where N is the total number of individuals. 
13 

14 

1 5 Spatial distributions of trees 
16 

* 

17 The unit-size m u -m u method and the p-index with successive changes of quadrat sizes (Iwao, 

18 1972) were used to analyze the spatial distribution of trees. Mean density m u is defined as: 

?„ 

19 m =J± — (10) 

20 where n . is the number of individuals in the jth quadrat and q u is the total number of quadrats 

21 when the quadrat size is u. The q u -values were 120, 60, 30, 15, 7, 4 and 2 respectively for the 
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1 w-values of 1 (2.5 m x 2.5 m), 2, 4, 8, 16, 30 and 60 in the Okinawa forest, and 64, 32, 16, 8, 4 

2 and 2 respectively for the w-values of 1 (2.5 m x 2.5 m), 2, 4, 8, 16 and 32 in the Ishigaki 

* 

3 forest. On the other hand, mean crowding m u is defined by Lloyd (1967) as: 



19 



* 



4 m u =^- q (11) 

2>, 

5 If the basic component of the spatial distribution is a single individual tree, individual 

6 trees are considered to be randomly distributed when m u = m u , aggregately distributed when 

7 m u > m u and uniformly distributed when m u < m u for any quadrat size. In order to provide 

8 knowledge on the distribution pattern of clumps, Iwao (1972) proposed the p-index. 

9 (12, 



* 

10 where for the smallest quadrat size (u = 1), p l = m l I m r When the values of p u are plotted 

1 1 against the quadrat sizes, a peak of the curve may suggest the clump area. 
12 

13 

14 Overlapping in spatial distributions of trees between layers 
15 

16 The a> -index, which is proposed by Iwao (1977) for analyzing spatial association between 

17 species, was applied to measure the degree of overlapping in spatial distributions of trees 

1 8 among layers with successive changes of quadrat size u. 



20 ^ = YV^ ir - r ^ 

ZU i_ /(ind) 



(13) 

10 
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1 



2 - Z. ^ (ind) 

( " ,= r„ c^W 

3 



4 where 7 and y (ind) are respectively given in the forms: 



5 A, 

6 



— — 2 



7 , Z"vT" 

V ;'=i y=i 
8 

1 7=1 7=1 



10 ,(md) a is 2 ^ ; 

11 h -' M 

12 Here, 2« is the total number of quadrats taken as 120, 60, 30, 15 and 7 respectively for the u- 

13 values of 1 (2.5 m x 2.5 m), 2, 4, 8 and 16 in the Okinawa forest, and 64, 32, 16, 8, 4 and 2 

14 respectively for the M-values of 1 (2.5 m x 2.5 m), 2, 4, 8 and 16 in the Ishigaki forest, riAj and 

1 5 ubj are the number of individuals of the j'th quadrat respectively belonging to Layer A and 

16 Layer B. The value of a> changes from the maximum of +1.0 for complete overlapping, 

17 through 0.0 for independent occurrence, to the minimum of -1 .0 for complete exclusion. 
18 

19 

20 Dendrogram 
21 

22 The dendrograms for analyzing the degrees of floristic similarity among layers were 

23 constructed following Mountford's (1962) method, or unweighted pair-group method using 

24 arithmetic averages (Sneath & Sokal, 1973). 
25 
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1 

2 Regression analysis 

3 

4 The coefficients for nonlinear equations were determined with statistical analysis software 

5 (KaleidaGraph V. 4.0, Synergy Software, USA). On the other hand, the coefficients for linear 

6 and curvilinear equations were determined with the least-squares method. 
7 

8 

9 RESULTS 

10 

1 1 Architectural stratification 
12 

13 The M-w diagram in the Okinawa forest is illustrated in Fig. la. As is clear from Fig. lb, the 

14 M-w diagram shows four phases, the first and fourth of which phases possess a property of 

15 Eq. (1), i.e. Type I of the C-D curve tribe that contains eight types (Shinozaki & Kira 1961), 

16 whereas the second and third phases possess a property of Eq. (2), i.e. a power function. The b 

17 values in Eq. (2) were significantly different between the second and the third phases (t = 212, 

18 P = 0). As a result, it was confirmed that the forest consists of four architectural layers. 

19 Figure lc shows that the extrema, represented by arrows, emerged apparently in the S 2 - In w 

20 relationship. Their tree weights at boundaries between layers were estimated to be 508, 3.18 

21 and 0.00953 cm 2 m. 
22 

23 On the other hand, the M-w in the Ishigaki forest is illustrated in Fig. 2a. As is clear from Fig. 

24 2b, the M-w diagram shows five phases, all of which have the property of Eq. (1). As a result, 

25 it was confirmed that the forest consists of five architectural layers. Figure 2c shows that the 
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1 extrema, represented by arrows, apparently emerged in the 5 2 - In w relationship. Their tree 

2 weights at the boundaries between layers were estimated to be 706.3, 18.74, 0.6871 and 

3 0.01236 cm 2 m. 
4 

5 Figure 3 shows the relationships between tree height H and weight w respectively for the 

6 Okinawa and Ishigaki forests, which were respectively formulated (cf. Kira & Ogawa, 1971) 

7 as follows: 

8 — = ^ + — ( 14 ) 

H 1.37w 0 ' 361 21.0 

9 77 = 1.06w 0368 + 2L0 (15) 

10 The estimated maximum tree height was 21.0 m, which was the same in both of the forests. 

1 1 The heights of the boundaries were determined as 8.0, 1.9 and 0.25 m in the Okinawa forest 

12 and 7.6, 2.7, 0.87 and 0.21 m in the Ishigaki forest by substituting the tree weights at 

13 boundaries obtained above for w respectively in Eqs. 14 and 15. Therefore, in the Okinawa 

14 forest, the height range was 8.0 m < H < 16.5 m for the top layer, 1 .9 m < H < 8.0 m for the 

15 second layer, 0.25 m < H < 1.9 m for the third layer and 0. 10 m < H < 0.25 m for the bottom 

16 layer. In the Ishigaki forest, the height range was 7.6 < H < 14. 1 m for the top layer, 2.7 < H < 

17 7.6 m for the second layer, 0.87 < H < 2.7 m for the third layer, 0.21 < H < 0.87 m for the 

1 8 fourth layer and 0.0 < H < 0.2 1 m for the bottom layer. 
19 

20 

2 1 Species dominance 

22 

23 Table 1 lists importance value IV of five woody species for each layer in each of the forests in 

24 order of species rank, which was determined from IV in the total stand. In the subtropical 
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1 evergreen broadleaf forest in Okinawa Island, a total of 26 families, 43 genera, 60 species and 

2 4684 woody individuals were recorded. The most species rich family was Rubiaceae, which 

3 contained 12 species. Symplocos, Lasianthus and Ilex were the species rich genera, each of 

4 which contained five species. Out of the 60 species, only three species (5%) consisted of a 

5 single individual. Castanopsis sieboldii (Mak.) Hatusima appeared in all layers with the 

6 highest importance value, especially with a tremendously high value of 44% in the top layer. 

7 This phenomena indicate that C. sieboldii is the most dominant and climax species. Schima 

8 wallichii (DC.) Korth. was the second dominant species in terms of IV of 9.30% in the total 

9 stand. A quite high TV of 22.8% of S. wallichi in the top layer compared with the very low IV 

10 (ranging from 0.62 to 1.69%) in the lower three layers indicate the heliophilic nature of the 

1 1 species. 
12 

13 In the subtropical evergreen broadleaf forest in Ishigaki Island, a total of 33 families, 52 

14 genera, 77 species and 4157 woody individuals were recorded. The most species-rich family 

15 was Aquifoliaceae, which contained six species. Aquifoliaceae consisted of Ilex alone, which 

16 is the most species-rich genus. Out of the 77 species, only four species (5%) consisted of a 

17 single individual. Ardisia quinquegona Blume was the most dominant species in terms of the 

18 highest importance value in the total stand and in all layers, except the top layer (Table 1). 

19 This species may be a subcanopy species, i.e. not a climax species, because of its 

20 disappearance in the top layer. Castanopsis sieblodii (Mak.) Hatusima was a low rank species 

21 in terms of the low importance value in the total stand, though it is a climax species in the 

22 Ishigaki subtropical evergreen broadleaf forest (Foster et al., 1960). However, this species 

23 appeared in all layers with a considerably higher importance value, especially the second 

24 highest importance value in the top layers (Data not shown). 
25 
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1 

2 Species-area relationship 

3 

4 As shown in Fig. 4, the expected number of species increased, and then tended to be saturated 

5 with increasing number of quadrats for the Okinawa forest, whereas it increased with 

6 increasing number of quadrats for the Ishigaki forest. The relationships of the expected 

7 number of species S g to the number of quadrats q in each layer and the total stand were well 

8 approximated by the following equation (Ogawa, 1980; cf. Hagihara, 1995): 

9 06) 

10 where c and d are coefficients, and 5 max is the expected maximum number of species. 
11 

12 In the Okinawa forest, the expected maximum number of species was estimated to be 62 in 

13 the total stand. The S max increased from 26 in the top layer, through 46 in the second layer 

14 and 61 in the third layer, to 65 in the bottom layer. Therefore, the bottom layer contained the 

15 highest potential number of species. This result was different from the Ishigaki forest, where 

16 the expected maximum number of species was estimated to be 92 in the total stand. The 

17 fourth layer contained the highest potential number of species as the S max increased from 44 in 

1 8 the top layer, through 55 in the second layer and 77 in the third layer, to 90 in the fourth layer, 

19 but decreased to 88 in the bottom layer. 
20 

21 

22 Floristic similarity 

23 

24 The floristic similarities among layers of each forest were classified using dendrograms of 
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1 similarity index C n , as shown in Fig. 5. In the Okinawa forest, the strongest similarity in 

2 floristic composition was marked between the second and the third layers with a C n -value of 

3 0.80. The next was between the second-third and the bottom layers with a C n -value of 0.64. 

4 The lowest C n -value of 0.48 was between the top and the lower three layers. On the other 

5 hand, in the Ishigaki forest, the third and the bottom layers showed the highest similarity in 

6 floristic composition with a C n -value of 0.94. The second and third highest similarities were 

7 respectively between the third-bottom and the fourth layers with a C n -value of 0.69, and 

8 between the third-bottom-fourth and the second layers with a C n -value of 0.36. The lowest 

9 C n -value of 0. 15 was between the top and the lower four layers. 
10 

11 

1 2 Woody species diversity in the stratified forest stands 
13 

14 In the Okinawa forest, the values of //'and /'were respectively 2.75 bit and 0.66 for the top 

15 layer, 4.37 bit and 0.81 for the second layer, 4.73 bit and 0.80 for the third layer, 4.33 bit and 

16 0.73 for the bottom layer. In the Ishigaki forest, on the other hand, the values of H' and J' 

17 were respectively 3.84 bit and 0.90 for the top layer, 4.49 bit and 0.87 for the second layer, 

18 4.29 bit and 0.71 for the third layer, 4.21 bit and 0.70 for the third layer, 3.73 bit and 0.65 for 

19 the bottom layer (Table 2). These results indicate that the values of diversity indices //' and 

20 7'tended to increase from the top layer downward except for the bottom layer in the Okinawa 

2 1 forest, whereas in the Ishigaki forest they tended to increase from the bottom layer upward 

22 except for the H ' -value of the top layer. 
23 
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1 When species diversity compared among layers, the highest value of H' was in the third 

2 layer in the Okinawa forest. This is because of the second highest value of /' and the highest 

3 species richness (56 species) in the third layer. However, in the Ishigaki forest, the highest 

4 value of H ' occurred in the second layer. This is mainly due to the second highest value of 

5 /' ^ though species richness (35 species) was the fourth highest in the second layer. On the 

6 other hand, the lowest value of H ' was in the top layer in the Okinawa forest. This is ascribed 

7 to the lowest /'-value and species richness (18 species) in the top layer. The lowest value of 

8 //'was in the bottom layer. This is ascribed to the lowest value of /' , though species 

9 richness (52 species) was the third highest. In addition, in the Okinawa forest, the values of 

10 //'in the second layer and in the bottom layer were nearly the same, while the value of 

1 1 J 'was higher in the second layer than in the bottom layer. This is because an increase of 

12 species richness from the second layer (41 species) to the bottom layer (53 species) could 

13 compensate for a decrease of evenness /' from the second layer to the bottom layer. On the 

14 other hand, in the Ishigaki forest, the values of //'in the top layer and in the bottom layer 

15 were nearly the same, while /'was the highest in the top layer, but was the lowest in the 

16 bottom layer. This is because an increase of evenness /' from the bottom layer to the top 

17 layer could compensate for a decrease of species richness from the bottom layer (52 species) 

18 to the top layer ( 1 9 species ) . 
19 

20 

21 Spatial distributions of trees 

22 

23 The spatial distribution patterns of trees for each layer are shown in Fig. 6 based on the unit- 

* 

24 size m u -m u relation and the p-index against quadrat size u. The spatial distribution of the 

25 basic component for all layers, except the top layer in the Okinawa forest and the bottom 

17 
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1 layer in the Ishigaki forest, was a single individual and its distribution was random, because 

* 

2 the m u -m u relation line mimicked the Poisson line at a significance level of 0.01. However, 

3 in the top layer of the Okinawa forest, the spatial distribution of trees showed a special trend. 

4 The p-index showed three peaks respectively at quadrat size 1 , 4 and 30. The first peak at 

5 quadrat size 1 may be occasional appearance due to the underdispersion of individual trees, 

6 and the second and third peaks may perhaps be related to topography. As a result, the top 

7 layer probably consists of a double clump, the small and large clumps of which were 

8 respectively 25 m 2 and 187.5 m 2 . In light of mean area occupied by individuals of the top 

9 layer trees in the Okinawa forest (9.1 + 0.6 (SE) m 2 ), the small clump may include three 

10 individuals and the large clump may include around eight small clumps. On the other hand, in 

1 1 the bottom layer of the Ishigaki forest, the spatial distribution of trees was aggregated (P = 1.1 

12 x 10" 4 ), though its basic component was a single individual. 
13 

14 

1 5 Overlapping in spatial distributions of trees among layers 
16 

17 The degree of overlapping co with successive changes of quadrat sizes in spatial 

18 distributions of trees among layers combined successively from the top layer downward is 

19 shown in Fig. 7. In the Okinawa forest, the spatial distributions of trees were overlapped 

20 between the top and the second layers, between the top-second and the third layer, and also 

21 between the top-second-third and the bottom layer. Similarly, in the Ishigaki forest, the spatial 

22 distributions of trees were nearly overlapped between the top-second and the third layers, 

23 between the top-second-third and the fourth layers, and also between the top-second-third- 

24 fourth and the bottom layers. Only a different degree of overlapping occurred between the top 

25 and the second layers, where the spatial distributions tended to be independent with each other. 
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1 These results may indicate that trees in the upper two layers in the Ishigaki forest can catch 

2 sufficient light, while light can not penetrate easily to the lower three layers in both of the 

3 forests. As a result, almost species in the lower layers might be shade-tolerant in both of the 

4 forests. 
5 

6 

7 Mean tree weight and density among layers 

8 

9 As shown in Fig. 8, mean tree weight w, of the ith layer decreased from the top (i= 1) toward 

10 the bottom layer (i = 4), whereas the opposite trend was observed for the tree density p t of the 

1 1 ith layer. This trend was successfully expressed for the forest in Okinawa Island in the form 

12 (Feroz et al., 2006): 



13 Wi =Kp- 



(17) 



PoJ 

14 In the Okinawa forest, the values of coefficients K , a, and p 0 were estimated as 2.23 x 10 8 

15 cm 2 m ha"", 1.49 and 24575 ha -1 . The relationship between mean tree weight w, and tree 

16 density p i in the Ishigaki forest was well described by Eq. 15 with the same values of the 

17 coefficients obtained in the Okinawa forest (R 2 = 0.98). 
18 

19 

20 DISCUSSION 

21 

22 All species in both of the forests in Okinawa Island and Ishigaki Island are evergreen 

23 broadleaf. The number of species in the Okinawa forest was lower than that in the Ishigaki 

24 forest, though the study area in Okinawa Island was approximately twice as large as that in 

19 
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1 Ishigaki Island. Similarly, the potential number of species in the Okinawa forest (62) was 

2 lower than that in the Ishigaki forest (92). 
3 

4 In both of the forests, Castanopsis sieboldii (Mak.) Hatusima, Daphniphyllum teijmannii Zoll. 

5 ex Kurz, Neolitsea aciculata (Bl.) Koidz. and Dendropanax trifidus (Thunb.) Makino were 

6 typically facultative shade species (light-tolerant under high light conditions and shade - 

7 tolerant under low light conditions) as they appeared in all layers, especially Castanopsis 

8 sieboldii was the most dominant species in the Okinawa forest with its highest importance 

9 value in each layer (Table 1). Although Ardisia quinquegona Blume was the most dominant 

10 species in the Ishigaki forest, it was not typically a facultative shade species. This is because 

1 1 this species did not appear in the top layer, i.e. mainly shade-tolerant, growing from the 

12 bottom layer to the subcanopy layer (or second layer). Ardisia quinquegona also did not 

13 appear in the top layer in the Okinawa forest, though it was one of the dominant species. 

14 Schima wallichii (DC.) Korth. in the Okinawa forest was typically a sun species as it appeared 

15 in the top layer with a tremendously high importance value, while it had a very low 

16 importance value in the lower three layers. In the Ishigaki forest, however, this species only 

17 appeared in the lower three layers, indicating a shade-tolerant nature which is different from 

18 the Okinawa forest. 
19 

20 The similarity in floristic composition between the top layer and the lower three layers in the 

21 Okinawa forest was weak, i.e. approximately one-third of the species from the lower three 

22 layers may be able to grow into the top layer. In the Ishigaki forest, however, the floristic 

23 composition between the top layer and the lower four layers was almost exclusive, i.e. 

24 approximately one-sixth of the species from the lower four layers may be able to grow into 

25 the top layer. This is because the floristic composition was more similar among layers in the 
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1 Okinawa forest than in the Ishigaki forest (Fig. 5). Moreover, a low similarity in floristic 

2 composition was observed between these two forests, which is concluded from the low value 

3 of similarity index ( C n = 0.33), though approximately half of the total number of species were 

4 common. This is mainly due to the different number of individuals belonging to the same 

5 species between these two forests. For example, the number of individuals of Castanopsis 

6 sieboldii (Mak.) Hatusima in the Okinawa forest (10480 ha" 1 ) was not the same as in the 

7 Ishigaki forest (2600 ha" 1 ). 
8 

9 It is known that the diversity of a community depends on two things: species richness and the 

10 evenness with which the individuals are apportioned among them (Pielou, 1975). As the 

1 1 lower layers contained many species relative to their smaller height ranges (Table 2), 

12 obviously these layers support high species richness of the forests. For example, 88% of the 

13 total species with 29% of the total individuals in the Okinawa forest and 67% of the total 

14 species with 30% of the total individuals in the Ishigaki forest packed within such thin bottom 

1 5 layers of 1 5 cm and 2 1 cm deep, respectively. 
16 

17 In the Okinawa forest, the value of H' for small-sized trees having H > 0.10 m was quite 

18 high as compared to that of H' for large-sized trees having DBH > 4.5 cm (Table 3). This is 

19 mainly caused by a large number of species for small-sized trees, though higher J' -value has 

20 a small influence on the high value of H' . As shown in Fig. 9, the trend of increasing 

21 diversity with successively decreasing height of layers from the top downward represents that 

22 high woody species diversity depended on small-sized trees. Since small-sized trees provide a 

23 natural habitat for animals living on the forest floor, conservation of small-sized trees in the 

24 lower layers is indispensable to sound maintenance of Okinawan evergreen broadleaf forests 

25 in Okinawa Island. 
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1 

2 On the other hand, in the Ishigaki forest, the value of H' for large-sized trees having DBH > 

3 4.5 cm was quite high as compared to that of H ' for the total stand (Table 3). This is mainly 

4 caused by a high evenness for large-sized trees, though smaller number of species has an 

5 influence on the low value of H' . The diversity tends to increase up to the second layer and 

6 then decrease downward (Fig. 9), which was different from the Okinawa forest. These results 

7 may represent that the large-sized trees have an important role in maintaining high woody 

8 species diversity. Thus, it may be thought that the decreasing trend of H' from the second 

9 layer downward is a characteristic of the subtropical evergreen broadleaf forest in Ishigaki 
10 Island. 

11 

12 The values of the diversity indices (//'and /') for the total stand in the Okinawa forest were 

13 higher than those in the Ishigaki forest (Table 3). This is because of a quite high evenness in 

14 the Okinawa forest compared to that in the Ishigaki forest, though species richness was lower 

15 in the Okinawa forest than in the Ishigaki forest. However, woody species diversity for large - 

16 sized trees is lower in the Okinawa forest than in the Ishigaki forest. This is caused by a 

17 higher species richness and evenness in the Okinawa forest than in the Ishigaki forest. 

18 Therefore, the statement of increasing multi-layering structure and woody species diversity 

19 along a latitudinal thermal gradient from higher latitudes to the tropics (Hozumi, 1975; 

20 Yamakura, 1987; Kira, 1991; Ohsawa, 1995; Kimmins, 2004) supported the present study for 

21 large-sized trees, but it did not support for total trees. 
22 

23 In Table 3, tree density of 62452 ha -1 for trees having H > 0.10 m in the Okinawa forest was 

24 quite low compared to that of 103925 ha -1 for total trees in the Ishigaki forest, whereas basal 

25 area of 50.3 m 2 ha" 1 in the Okinawa forest was almost the same as that of 49.4 m 2 ha" 1 in the 
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1 Ishigaki forest. However, tree density and basal area for trees having DBH > 4.5 cm in the 

2 Okinawa forest were respectively 3080 ha~'and 43.2 m 2 ha -1 , whose values are comparable 

3 with 2775 ha _1 and 42.7 m 2 ha _1 in the Ishigaki forest. Therefore, it follows that the Ishigaki 

4 forest is densely populated with small-sized trees compared to Okinawa forest. However, the 

5 biomasses of these forest stands are possibly similar, because the basal areas are very close to 

6 each other for trees having DBH > 4.5 cm. 
7 

8 The trend that mean tree weight decreased from the top toward the bottom layer, whereas tree 

9 density increased from the top downward (Fig. 8) was successfully formulated with Eq. 15. 

10 This trend resembles the mean plant weight-density trajectory of self -thinning even-aged 

1 1 plant populations (Hagihara, 2000), which start growing from initial plant densities lower than 

12 the initial plant density of the population obeying the -3/2 power law of self-thinning (Yoda et 

13 al., 1963). The relationship of mean tree height to tree density for the upper two layers in the 

14 Okinawa forest supported Yamakura's quasi -1/2 power law of tree height, because the a- 

15 value in Eq. 15 was close to 3/2 (Feroz et al., 2006). Our formulated equation is more general 

16 than Yamakura's quasi -1/2 power law of tree height, because the equation describes lower 

17 layers as well. 
18 

19 
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Table 1. Five dominant species in the Okinawa forest (a) and the Ishigaki forest (b) in order of species rank determined by importance value TV 
in the total stand. 



(a) 



Species 


Name of species 






IV % 






rank 




Top layer 


Second layer Third layer 


Bottom layer 


Total stand 


1 


Castanopsis sieboldii (Mak.) Hatusima 


43.60 


12.63 


26.12 


21.74 


21.30 


2 


Schima wallichii (DC.) Korth. 


22.34 


1.78 


0.62 


1.69 


9.30 


3 


Elaeocarpus japonicus Sieb. & Zucc. 


10.16 


7.56 


1.73 


2.78 


4.61 


4 


Myrsine seguinii Lev'l. 




4.14 


5.46 


7.01 


3.92 


5 


Ardisia quinquegona Bl. 




4.97 


6.88 


3.72 


3.74 


(b) 


Species 


Name of species 






IV % 






rank 




Top layer 


Second 


Third Fourth Bottom 


Total 








layer 


layer layer 


layer 


stand 


1 


Ardisia quinquegona Bl. 




9.23 


24.66 10.52 


24.71 


15.39 














2 








8.53 19.10 


9.45 


7.93 




Lasianthus plagiophyllus Hance 






3 








9.34 6.89 


6.88 


6.54 




Psychotria manillensis Bartl. 






4 






0.84 


0.79 10.98 


3.18 


6.19 




Sarcandra glabra (Thunb.) Nakai 




5 


Camellia sasanqua Thunb. 




7.55 


1.65 6.82 


6.73 


5.54 







28 



Downloaded from http://biorxiv.org/on September 18, 2014 



1 Table 2. Diversity indices among layers in the Okinawa forest (a) and the Ishigaki forest (b). 
2 

3 (a) 



Layer 


Height range (m) 


No. of trees No. of species 


H' (bit) 










per plot area 






Top 


8.0 <#< 16.5 


139 


18 


2.75 


0.66 


Second 


1.9<#<8.0 


953 


42 


4.37 


0.81 


Third 


0.25 <H< 1.9 


2226 


56 


4.73 


0.80 


Bottom 


0.10<//<0.25 


1366 


51 


4.33 


0.73 


(b) 


Top 


7.6 <#< 14.1 


42 


19 


3.84 


0.90 


Second 


2.7 <H<1.6 


130 


35 


4.49 


0.87 


Third 


0.87 <H< 2.7 


915 


56 


4.29 


0.71 


Fourth 


0.21 <H< 0.87 


1806 


65 


4.21 


0.70 


Bottom 


0.0 <//< 0.21 


1264 


52 


3.73 


0.65 



6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
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Table 3. Comparison in species diversity and tree density between the Okinawa and Ishigaki forests. 



Study site 


Sample 


Tree size 


Basal area No. of 


S 


H' (bit) 


J' 


Source 




area (ha) 




(m 2 ha _1 ) 


trees (ha" 










Okinawa forest 


0.075 


#>0.10m 
DBH>4.5 cm 


50.3 
43.2 


62452 
3080 


60 
32 


4.83 
3.80 


0.82 
0.77 


Feroz et al. 2006 


Ishigaki forest 


0.04 


#>0.0m 


49.4 


103925 


77 


4.36 


0.69 








DBH>4.5 cm 


42.7 


2775 


36 


4.61 


0.89 





S, the number of species per sample area. 
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1 FIGURE LEGENDS 

2 

3 Fig. 1. Relationships of mean tree weight M (a), the first derivative Si (b) and the second 

4 derivative 52 (c) to tree weight w on logarithmic coordinates in the Okinawa forest. In the 

5 M-w diagram (a), the regression curves of the top and bottom layers are given by equation 

6 (1) (R 2 = 0.99) for the top layer and (R 2 = 0.97) for the bottom layer. The regression 

7 curves of the second and third layers are given by equation (2) (R 2 = 0.91) for the second 

8 layer and (R 2 = 0.94) for the third layer. The Si and S2 are respectively defined as 

9 equations (3) and (4). The moving values of Si and S2 are plotted against their 

10 corresponding geometric means of w. 

1 1 Fig. 2. Relationships of (a) mean tree weight M, (b) the first derivative Si and (c) the second 

12 derivative S2 to tree weight w on logarithmic coordinates in the Ishigaki forest. In the M-w 

13 diagram (a), the regression curves for all layers are given by Eq. (1) (R 2 = 0.98) for the top 

14 layer, (R 2 = 0.92) for the second layer, (R 2 = 0.90) for the third layer, (R 2 = 0.88) for the 

15 fourth layer and (R 2 = 0.96) for the bottom layer. The application of Si and S2is the same 

16 as in Fig. 1 . 

17 Fig. 3. Relationships between tree height H and weight w. In the Okinawa forest (a), the 

18 regression curve is given by Eq. (14) {R 2 = 0.97). Filled circles, top layer (8.0 m < H < 

19 16.5 m); open circles, second layer (1.9 m < H< 8.0 m); Filled squares, third layer (0.25 

20 m < H < 1.9 m); open squares, bottom layer (0. 10 m <H< 0.25 m). In the Ishigaki forest 

21 (b), the regression curve is given by Eq. (15) (R 2 = 0.96). Filled circles, top layer (7.6 m < 

22 H < 14. 1 m); open circles, second layer (2.7 m < H < 7.6 m); Filled squares, third layer 

23 (0.87 m <H< 2.7 m); open diamond, fourth layer (0.21 m < H < 0.87 m); open squares, 

24 bottom layer (0.0 m < H < 0.21 m). 

25 Fig. 4. Species-area curves (area of one quadrat: 2.5 m X 2.5 m). 
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1 Filled diamond, total stand and other symbols are the same as in Fig. 3. In the Okinawa 

2 forest (a), the curves are given by equation (16), where Smax = 62.4 (R 2 = 1.0) for the total 

3 stand, S max = 32.8 (R 2 = 1 .0) for the top layer, S max = 47.0 (R 2 = 1 .0) for the second layer, 

4 Smax = 60.30 (R 2 = 1.0) for the third layer, S max = 61.1 (R 2 = 1.0) for the bottom layer. In 

5 the Ishigaki forest (b), the curves are given by the same equation, where S max = 92.6 (R 2 = 

6 1 .0) for the total stand, S max = 43.7 (R 2 = 1 .0) for the top layer, S max = 54.7 (R 2 = 1 .0) for 

7 the second layer, 5 max = 7 1 .5 (R 2 = 1 .0) for the third layer, S max = 90.2 (R 2 = 1 .0) for the 

8 fourth layer, 5 max = 77.9 (R 2 = 1 .0) for the bottom layer. 

9 Fig. 5. Dendrograms of the degree of similarity C n in floristic composition among layers in 

10 the Okinawa forest (a) and in the Ishigaki forest (b). 

* 

1 1 Fig. 6. A schematic representation showing the relationships between mean crowding m u and 

12 mean density m u , and p -index with successive changes of quadrat size u in the Okinawa 

13 forest (a) and in the Ishigaki forest (b). 

14 Fig. 7. Degree of overlapping a> with successive changes of quadrat size u in spatial 

15 distributions of trees in the Okinawa forest (a) and in the Ishigaki forest (b). The smallest 

16 quadrat size (u = 1) is 2.5 m X 2.5 m. 

17 Fig. 8. Relationship between mean tree weight w, and tree density p i among the layers. The 

18 curves are given by Eq. (17) in the Okinawa and Ishigaki forests (R 2 = 0.98). Filled circle, 

19 Okinawa forest; open circle, Ishigaki forest. 

20 Fig. 9. Trends of H' along the layers combined from the top downward in the Okinawa forest 

21 (a) and in the Ishigaki forest (b). Symbols are the same as in Fig. 8. 
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Figure 5 
Feroz et al. 
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Figure 6 
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Figure 9 
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